home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
IRIX 6.2 Development Libraries
/
SGI IRIX 6.2 Development Libraries.iso
/
dist
/
complib.idb
/
usr
/
share
/
catman
/
p_man
/
cat3
/
complib
/
DSISL.z
/
DSISL
Wrap
Text File
|
1996-03-14
|
2KB
|
67 lines
DDDDSSSSIIIISSSSLLLL((((3333FFFF)))) DDDDSSSSIIIISSSSLLLL((((3333FFFF))))
NNNNAAAAMMMMEEEE
DSISL - DSISL solves the double precision symmetric system A * X = B
using the factors computed by DSIFA.
SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
SUBROUTINE DSISL(A,LDA,N,KPVT,B)
DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
On Entry
AAAA DOUBLE PRECISION(LDA,N)
the output from DSIFA.
LLLLDDDDAAAA INTEGER
the leading dimension of the array A .
NNNN INTEGER
the order of the matrix A .
KKKKPPPPVVVVTTTT INTEGER(N)
the pivot vector from DSIFA.
BBBB DOUBLE PRECISION(N)
the right hand side vector. On Return
BBBB the solution vector X . Error Condition
AAAA division by zero may occur if DSICO has set RCOND .EQ. 0.0 or DSIFA
has set INFO .NE. 0 . To compute INVERSE(A) * C where C is a matrix
with P columns
CALL DSIFA(A,LDA,N,KPVT,INFO)
IF (INFO .NE. 0) GO TO ...
DO 10 J = 1, P
CALL DSISL(A,LDA,N,KPVT,C(1,J)) 10 CONTINUE LINPACK. This version
dated 08/14/78 . James Bunch, Univ. Calif. San Diego, Argonne Nat. Lab.
Subroutines and Functions BLAS DAXPY,DDOT Fortran IABS
PPPPaaaaggggeeee 1111